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As an application of perfect lattice perturbation theory, we construct 
an 0(A) perfect lattice action for the anharmonic oscillator analytically in 
momentum space. In coordinate space we obtain a set of 2-spin and 4-spin 
couplings oc A, which we evaluate for various masses. These couplings never 
involve variables separated by more than two lattice spacings. 

The 0(A) perfect action is simulated and compared to the standard ac- 
tion. We discuss the improvement for the first two energy gaps AEi, AE2 
and for the scaling quantity AE2/AE1 in different regimes of the interaction 
parameter, and of the correlation length. 
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1 Introduction 



The only non-perturbative access to complicated 4d quantum field theories, 
such as QCD, which proved successful, are Monte Carlo simulations on the 
lattice. They necessarily take place at a finite lattice spacing a and in a 
finite size L. In order to reveal information about continuum physics in an 
infinite volume, we have to require a << ^ << L, where ^ is the correlation 
length. In particular the finiteness of ^/a causes serious systematic errors 
in practical simulations. It is now very fashionable to fight such artifacts 
by using "improved lattice actions" These are discretizations of the 
continuum action, which are supposed to display the correct scaling behavior 
down to a much shorter correlation length in lattice units, than it is the case 
for the standard lattice action. 

In the literature, there are mainly two strategies to construct improved 
lattice actions, in particular for QCD. The first one is called Symanzik's pro- 
gram 0. One tries to eliminate the lattice spacing artifacts order by order 
in a - similar to the Runge and Kutta procedure for the numerical solution 
of ordinary differential equations. This is achieved by adding irrelevant oper- 
ators. On the classical level, and in the framework of on-shell improvement 
[^], the standard Wilson action for QCD could be improved to 0{a) ana- 
lytically by adding the so-called clover term 0. On the quantum level, the 
coefficient of this term gets renormalized, and the quantum correction was 
first estimated numerically by a mean field approach The complete 0{a) 
improvement was finally determined by the ALPHA collaboration based on 
extensive simulations However, it seems hardly feasible to carry on this 
program beyond 0(a). 

The alternative method uses renormalization group concepts to construct 
quasi-perfect actions. These are approximations to perfect actions, i.e. to 
actions which are completely free of cutoff artifacts 0]. As a fundamental 
difference from Symanzik's program, this method is non-perturbative with 
respect to a. As a first step, this program can be realized perturbatively (in 
the interaction), which yields analytic expressions for the perfect quark-gluon 
and 3-gluon vertex functions [|, |10|- Q Thus one eliminates all artifacts 
of 0(a") and 0{ga"'), such that the remaining artifacts are of 0{g'^a) and 
beyond {g is the gauge coupling). This is opposed to the action of Ref. 

^Perturbatively perfect actions have also be studied for the Schwinger model, jl^, [l2| . 



2 



which is free of artifacts in 0{g"a), but plagued for instance by systematic 
errors in O(a^). 

An extension of this program is the construction of "classically perfect 
actions" [|13|. This approach, which is designed particularly for asymptoti- 
cally free models, is non-perturbative also with respect to the coupling g. 
Using a multigrid procedure, one identifies the fixed point action of an renor- 
malization group transformation. This can be done solely by minimization 
- the functional integral reduces to a classical field theory problem - and 
the fixed point action then serves as an approximatively perfect ( "classically 
perfect" ) action at finite correlation length too. In a sequence of toy models, 
it turned out that classically perfect actions are excellent approximations to 
(quantum) perfect actions, in the sense that they drastically suppress lattice 
spacing artifacts. The improvement achieved in this way goes far beyond first 
order Symanzik improvement. This has been observed for the 2d 0(3) and 
CP(3) model |T^, the Schwinger model |]T5| and the Id XY model 



In principle that program can be extended also beyond classical perfection, 
if one performs e.g. one real space MCRG step at finite correlation length, 
starting from a classically perfect action. 

The construction, which is non-perturbative in a and in g, is presum- 
ably the climax of the improvement program. However, in perfect and also 
in classically perfect actions the couplings tend to involve infinite distances, 
and we can at best achieve locality in the sense of their exponential decay. 
For practical purposes a truncation is needed, which does some harm to the 
quality of the improvement. This is the main reason why the second, more 
sophisticated, improvement program could not be applied yet in a satisfac- 
tory way to QCD. 

Here we focus on the perturbatively perfect action. It has potential appli- 
cations with two respects: it can either be used directly, or as a starting point 
of the non-perturbative multigrid improvement [0. A direct application of 
a truncated perfect quark-gluon vertex function - together with truncated 
perfect free quarks - to heavy quarks is presently under investigation. Pre- 
liminary results for the charmonium spectrum are given in Ref. . 



The purpose of this paper is to test specifically such a direct application 
in a very simple situation. Our model is the Id model, or anharmonic 
oscillator. 

As a toy model, the anharmonic oscillator has a number of virtues: we 
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can achieve an excellent locality, such that our 0(A) perfect action does not 
need any truncation of the couplings. Thus the perturbative improvement 
can be tested separately, without admixture of truncation effects. Moreover, 
our construction is based on continuum perturbation theory, and there we 
do not encounter any divergent loop integrals, in contrast to field theory 
{d > 1). Finally, the reduction to quantum mechanics has the advantage 
that the couplings we identify do not get renormalizcd in the full theory. 

On the other hand, it is exceptionally difficult in our model to demon- 
strate an improvement compared to the standard action, because the latter 
is also very good in this case. For the harmonic oscillator it is even perfect 
itself, for small interaction - the regime of interest here - it is still very good, 
and even for moderate interactions it performs amazingly well. As a further 
problem we note that the performance of continuum perturbation theory, 
which our improvement is based on, is rather poor in this model. 

The advantages and disadvantages listed above are specific for the one 
dimensional case. 



2 The model in the continuum 

The observables we are going to consider can be evaluated directly in the con- 
tinuum to a fantastic accuracy. Our interest is of course not in their values, 
but solely in the comparison of lattice artifacts in different discretizations. 
We want to test the success of a specific improvement program for the lattice 
action. 

Nevertheless we have to start by recalling some properties of the conti- 
nuum system. To fix the (field theoretic) notation, we denote the Euclidean 
action as 

^Jdt [^m' + + Mtr] ■ (2.1) 

Throughout this paper we assume m, A > 0, hence we only study the "sym- 
metric phase" (as opposed to the double well). We consider the energy 

eigenvalues En, more precisely we are going to measure directly the energy 
gaps AEn = En — Eq. An additive constant all over the spectrum is out of 
control, and not much of interest either. The simplest scaling quantity is 

f|^A£,.«, (2.2) 



4 



where ^ is the correlation length. Moreover, there is the simple relation 



AE„K,A) ^^y/^>0- (2.3) 

This just follows from rescaling t ^ t/fi and momentum k —>■ nk in the 
Hamiltonian (noi rescaling all dimensional quantities). 

The quantity ( |2.2D is a scaling quantity in the strict sense, i.e. a dimen- 
sionless ratio of physical observables. On the other hand, quantities like ( |2.3| ) 
may also involve unphysical normalization factors. In a field theoretic lan- 
guage they correspond to the asymptotic scaling. Actually, improved actions 
are designed for an improvement of scaling, but the influence on asymptotic 
scaling is of interest too. It has been observed before for the Gross Neveu 
model ||2^ and for pure SU(3) gauge theory |2l| that "accidentally" the latter 



is also improved for (quasi-)perfect actions. 

For m > 0, the strength of the interaction depends on the dimensionless 
parameter 

A = A ' (2-4) 

which is obvious from eq. (|2.3| ). The energy eigenvalues can be expanded 
in A, but these expansions diverge at large orders (the coefficients oscillate 
and their absolute values grow faster than any polynomial). Q However, 
a truncated series is still useful at A << 1 (this situation is familiar from 
QED). The coefficients of these expansions have been derived many times 



in the literature, for instance in Ref. 22 



^(A) ^ i + 3A_i8p+ ™P-3825A^ 
m 8 

m ^ ' 4 8 64 ' 

^(X).2.3A-lfP.IflA3-155|5£x', (2.5, 



^This argument was made rigorous first by K. Symanzik (unpublished). The point is 
that the rescahng can be implemented unitarily jl^ . 

■^Note that the point A = is non analytic. For a discussion of large orders, see e.g. 
Rcf. Hi. 

''^We expect the same behavior also for the the couplings in the perturbatively perfect 
action. 



5 



Tables of explicit values at finite A are given for example in Refs. |T^, ^ 
However, in particular for small A they can easily be reproduced from an 
eigenvalue problem, as described for instance in Ref. ||2^ . 

Since our construction in the following sections is also perturbative in A, 
it is important to know how the perturbation series for the above quantities 
behave. Figs. |l] and ^ compare the exact function to the truncated expansion 
in first, second, third and fourth order. 
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Figure 1: The ratios AEi/m (on top) and AE2/m (below) as functions of X. 
The exact result is compared to the perturbation series truncated at various 
orders. 
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Figure 2: The ratio AE2/AE1 as a function of X. The exact result is com- 
pared to the perturbation series truncated at various orders. 

We see that the apphcabihty of truncated expansions is restricted to very 
small values of A. The exact range depends on the quantity considered; it 
agrees with the relative magnitude of the coefficients in the expansions ( p.5|) . 
This range gradually expands as we proceed to higher orders. 

3 The 0(A) perfect action in momentum space 

If a system is given by some lattice action, then its physical properties re- 
main unaltered under a block variable renormalization group transformation 
(RGT) For suitable RGT parameters and infinite correlation length, an 
infinite number of iterations may lead to a finite fixed point action (FPA). 
A FPA is an example for a perfect action, since it is invariant under an 
RGT, and hence insensitive to the lattice spacing. Perfect actions also exist 
at any finite correlation length 0]. They reveal exact continuum scaling at 
any lattice spacing. For free or perturbatively interacting fields, they can be 
computed analytically in momentum space. This calculation simplifies if we 
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send the blocking factor to infinity, which amounts to a technique that we 
call "blocking from the continuum" . Recently this method has been applied 
to the Schwinger model |]rT| and to QCD Here we want to apply it to 
construct a lattice action for the anharmonic oscillator, which is perfect to 
0(A). Our blocking uses the standard, piece-wise constant weight distribu- 
tion for the original variables, and a Gaussian transformation term. For an 



alternative ansatz, closer to the spirit of a decimation RGT, see Ref. |25 



Our perfect lattice action S[(j)] is determined by the functional integral 

R[4>,V] = :^ E - r'^'dt ^(t))' . (3.1) 

Here (p is the lattice field, x are the lattice sites and the continuum action 
s[ip] is given in eq. ( p.l| ). The RGT parameter a > is arbitrary; for any 
value of a the RGT keeps the partition function invariant, 

Z = Jd^ e-'^I'^l cx e-^t-^], (D0 = J] / (3-2) 

and with it all expectation values. The limit a corresponds to the well 
known "5 function RGT" . 

In momentum space, this expression can be written as 

e-^t"^] = J D^Da exp { - ^ dk x 

[- E - 2nl)[{k + 2-Klf + m^]^{k + 27r/) 

2 lex 

Ma{-k)[<p{k) -Y,ip{k + 27(l)U{k + 2nl)] 

-l-^a a{—k)a{k)^ | x 

A /■ ^ 

{l - J d^P ^{Pi)^{P2MP3MP4) SiJ^Pi) 

+0(A^)}, 

H(A;) = p A; = 2sin^, (3.3) 
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where we have introduced an auxihary lattice field a. 
We denote the free continuum propagator as 



and 

G{k) = Y./\{k + 2Til)Ii{k + 2'Klf + a (3.5) 

is the perfect free lattice propagator, as we will see. At m = this is the 
fixed point propagator, which has been calculated for free scalar theories first 
by Bell and Wilson, and it characterizes the FPA for the O(A^) model in the 
large N limit as well |^ . 



We now choose the special value a = (sinhm — m)/m^, which renders 
the free lattice action "ultralocal" I^S], i.e. it only couples nearest neighbor 
lattice variables, ^ 

sinhm -m^ 1 ^ = 2sinhl^. (3.6) 

Our first step is the substitution 

^{k + 2tiI) = ip{k + 2tiI) - ia{k)/\{k + 27rZ)n(A; + 27r/), (3.7) 

which allows us to integrate out the continuum variable (p. We omit the 
constant factor in the Gaussian integral and obtain 



J Daexp[-^ dk [ia{-k)(P{k) + ^a{-k)G{k)a{k)]} x 
+6^7^ fn E TT r dpicr{p,)A{pi + 27mi)U{p, + 2nn, 

Zni '- ■_i „ ^rj7 -^VT J— TT 



=1 niS^ 
X27r5(pi +]92)5ni -nz 



^This choice for the RGT parameter a also provides optimal locality in d = 4 p9|| . 

^This is an example of an uncontrolled additive constant in S'[0], which motivates the 
consideration of the energy ga'ps^ rather than the single eigenvalues. The same holds for 
the subsequent integration over cr, see below. 
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f n E cr{pi)A{p, + 27m,)n(p, + 2 

x27r5(^[p, + 27rn,])+0(A2)}. 



(3.8) 



In a sense, this computation goes beyond the perfect QCD vertex function 
of Ref. 1^, because it includes - for the first time in the construction of a 
perfect action - a loop calculation, i.e. a quantum correction. The continuum 
loop integral reads 



A(x) 



x=0 



2tt 



dk A{k) 



1 
2m 



(m > 0), 



(3.9) 



which has been inserted above. In field theory we would encounter diver- 
gences at this point, which could be regularized by some standard technique 
in the continuum. Here the expression is finite from the beginning, and we 
will see that even the limit m —>■ can safely been taken at the end, when 
we identify the couplings in the 0(A) perfect lattice action. 
After performing a second substitution. 



a{k) = a{k)+tG{k)-'(l){k), 
we can integrate a, and we arrive at a lattice action of the form 



(3.10) 



+X\A+^ r dk (P{-k)B{k)(j){k) 



27r 



+ 



{2'Ky- 



d^'pCipmp^mp^mpsmp,) +0i\'). (3.11) 



This confirms that G is the perfect free lattice propagator for the RGT chosen 
here. Since it differs from the standard propagator only by a constant factor 
and a transformation of the mass, we can also confirm the statement that 
for the harmonic oscillator the standard action is perfect already. [] This is 
very specific for the case d = 1 considered here. 



^In fact, any lattice action for the harmonic osciUator is perfect 
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The functions B{k) and C{p) = C{pi,p2-P?,,P4) represent additional 2- 
variable and 4- variable couplings, while A is a constant, which is not really 
of interest to 0{X). 

The Wick contraction of two isolated a variables yields the lattice loop 
integral 



= r G{k)-^/^{k^2'KlfTl{k^2T:lf 



1 



m ■ sinh m ■ im? 

1 7 11 2 ^/ 4N 

^ 7777 + 777777"^ + OinT), 

2m 30 630 ^ ^' 



2 cosh m + e '"(I + sinh m) sinh m 

m 



(3.12) 



which obeys 7(m) — 7(— m) — 1/m. 

If we just insert this everywhere, we obtain the following 0(A) terms. 



^0 

Bo{k) 



1 

2m 

^ - 2-f{m)]G{k)-^ Yl ^(^ + 27r/)^n(A; + 2T:lf 
3m2 



m 



\- - 27(^) 



X 



2(sinhm • m^y 
^ + m) + k^m^m + 2m^| 



where 

C{p) 



m 



sinhm\ 



[ n <^bO"' E + 27rni)n(pi + 27rni) 



i=l 
4 



X 



SiYlPi + '^T^rii]) . 



(3.13) 



Note, however, that in the contractions in the cr^ term, an expectation 
value {a{pi)a(pj)) only enforces Pi — —pj ('pairing' of the lattice momenta), 
while the related summation integers n^, nj remain independent (no 'pairing' 
of the continuum momenta). Therefore, such contractions yield an additional 
contribution S{k), {k g] — tt, tt]), which does not factorize, 

1 r 

dqG{q)-^[ Y A(5 + 27rni)n(5 + 27rn2) X 

ni,n2,n3,n4 



S{k) 



271 
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A(-g + 27rn2)n(-g + 27m2)A(A; + 27m3)n(A; + 27^3) 
A{-k + 2nn^)U{-k + 27in,)] 5^^„^^o (1 - (3.14) 

(The case ni+n2 = = 723 + 77,4 is excluded here, because it has been included 



before in eq. ( p.l3| ).) This modifies the constant and the bilinear term to 

A = Ao + ^ r dk G{k)-^S{k), 

271 J -IT 

B{k) = Bo{k)-6G{k)-^S{k), (3.15) 

while the vertex function C{p), given in eq. (|3.13| ), is not affected. It does not 
pick up any loop contributions, and its structure can easily be understood in 
the language of "building blocks" as introduced for the quark-gluon vertex 

i- 

The lattice spacing artifacts in the standard action can be classified in 
magnitudes as Xa^, Aa^, Aa^, . . . , A^a^, A^a"^, . . . , A^a^ .... In the action we 
have constructed now, all artifacts oc A are erased. The remaining artifacts 
can still be oc a^, but they are multiphed at least by A^. This is analogous 
to QCD with the perfect vertex function, where the gauge coupling g plays 
the role of A. There the artifacts of the Wilson action start even in 0{ga), 
and the perturbative perfection pushes the leading artifact to 0{g'^a). 



4 The perfect couplings in coordinate space 

The interaction terms involving B and C, which we derived in momentum 
space, turn into convolutions in coordinate space. Hence B{r) describes 
additional 2- variable couplings ( "2-spin couplings" in a solid state language) 
over a distance r G and C(ri, r2, r^, 0) = C(r) introduces 4-spin couplings, 

a;,rG^ xG^,rG^3 

We exploit lattice translational invariance, in the latter case by setting r4 
to the arbitrary value 0. C(ri,r2r3,0) is invariant under permutation of its 
components, and C(r) = C(— r) (but there is no invariance under sign fiip 
of just one or two components of r). 
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4.1 The 2-spin couplings 

The function B{k) is even, hence 



B{r) = Bo{r) + 5i(r) = ^ T dk [Bo{k) - 6G{k)-^S{k)] cos{kr). (4.1) 

The first term, Bo{r), can be computed analytically. It turns out to be the 
dominating contribution to B{r), 



Pi 
P2 



2m 



— 7(m) 



m 



X 



(sinhm ■ m^)^ 

2 (3 + m^)m + 3m'^ + 2m^ 
— (4 + rh?)m — 2m^, 



J}. 



m H — m . 
2 



(4.2) 



(The quantities m and fh are defined in eqs. ( p.6[ ) and (|3.13|) .) 

Expanding in small m, we recognize the finiteness of this expression in 
the limit m — > 0, 



5o(r) 



l(0) 



(0)r 



r,-2 



(0) 



6f 



77 

Too 

91 



419 

1400 
1637 



m2 + 0(m^), 
m2 + 0(m^). 



300 12600 

1 0(^4) 

600 5040 ^ ^ 



(4.3) 



The additional term, Bi (r) , has to be evaluated numerically. It is signifi- 
cantly suppressed, essentially because the case = 0, i = 1 ... 4 is excluded 
from the summation in eq. ( p.l4| ). It typically affects the bilinear couplings 
only in third digit. 

^The finiteness at m = ("quartic oscillator") is a very sensitive consistency test. We 
did not study that case - where A = oo and En = c„A^/^ - extensively, since our improved 
action is designed for small A. However, we observed that Eq oc X^^^ can be fitted better 
for the perturbatively perfect action than for the standard action. 
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It turns out that also in -Bi(r) the couphngs are restricted to distances 
< 2. Therefore, the entire bihnear term B{r) is given by bi = bf''^ + b^j^\ 
i = 0, 1, 2. These couphngs are shown as functions of the mass in Fig. and 
some precise values are given Table 0. For completeness we also include the 
constant A. 




0.2 0.4 0.6 0.8 1 



Figure 3: The perfect bilinear couplings oc A, as functions of the mass. 



In addition we have of course the bilinear couplings of the free theory 
(resp. harmonic oscillator), given in eqs. ( pl6|) and In coordinate 

space the free action reads 



S[<l>]x=o 



m 



2 sinh m ■ rh? 



</)a; (2 + m^)S^^y - - (j)y . (4.4) 



4.2 The 4-spin couplings 

We recall that we describe the 4-spin couplings by C(r), r = (ri, r2, r^) G Z^, 
r4 = 0, that we have permutation invariance among ri, . . . , r4 and invariance 
under r — > —f, hence 



C{r} 



1 

(2^ 



d^p C{p) cos(p ■ r) . 



(4.5) 
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III = 


//; = 0.2 


= 0.3 


A 


0.1634448589 


0.1627118789 


0.1615030089 


bo 


0.7648784176 


0.7560248574 


0.7415398190 


hi 


0.3034599619 


0.2995928547 


0.2932745430 


b2 


0.0112505782 


0.0110757933 


0.0107911914 


Co 


0.2242402815 


0.2202292735 


0.2137198561 


Ci 


0.0579948140 


0.0568473063 


0.0549896682 


C2 


0.0002458836 


0.0002396038 


0.0002295148 


Cn 


0.0386197986 


0.0378246484 


0.0365388182 


C12 


0.0008364070 


0.0008162208 


0.0007837264 


C22 


0.0000983088 


0.0000956674 


0.0000914309 


Cn2 


0.0042341197 


0.0041403776 


0.0039890713 





m = 0.4 


m — 0.5 


m — 1 


A 


0.1598369028 


0.1577387270 


0.14198250413 


bo 


0.7218143338 


0.6973665904 


0.52780256487 


bi 


0.2846878032 


0.2740741694 


0.20145092650 


b2 


0.0104063833 


0.0099339708 


0.00680856517 


Co 


0.2049621122 


0.1942810514 


0.12580113799 


Ci 


0.0524997264 


0.0494780115 


0.03056403944 


C2 


0.0002161447 


0.0002001630 


0.00010708640 


Cn 


0.0348181119 


0.0327344171 


0.01982579170 


C12 


0.0007405382 


0.0006887131 


0.00038135880 


C22 


0.0000858306 


0.0000791587 


0.00004091627 


Cn2 


0.0037871600 


0.0035435707 


0.00206185706 



Tabic 1: The 0(A) perfect couplings for masses . . . 1. 
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Note that all the singularities at = are removable, both, for finite 
and for vanishing mass. 

It turns out that again the couplings never involve any two spins separated 
by a distance larger than 2. A general argument for that is given in the 
following subsection. This means that there are just 7 independent 4-spin 
couplings. We denote them as 

Co = C(0); Ci = 0, 0); C2 = C(2, 0, 0) 

Cn = C(l,l,0); Ci2 = C(l,2,0); C22 = C(2, 2, 0); Cn2 = 1, 2). 

They all represent equivalence classes, which contain a total of 65 nontrivial 
couplings (keeping = fixed). Some exact values are given in Table [1], 
and their mass dependence is illustrated in Fig. ^. 

4.3 Locality 

Assume that we calculate the perfect action to 0(A") (n > 1). This involves 
a number of perturbative correction terms, which arise from the expectation 
values of the continuum field if to some power. The highest power is (v?^")- 
There, Wick contractions lead to several kinds of terms, which we classify by 
the power of the inverse free propagator G^^. The maximal power is 

For instance, in the bilinear term, which depends only on one momentum, 
the maximal factor is G{k)~'^^. Therefore the maximal power of k is A;^" oc 
(1 — cosA;)^", which can be decomposed into terms oc cos/c, . . . ,cos2nk. In 
coordinate space this yields couplings oc 6x,i + S^-i, . . . , [Sx,2n + Sx-2n]- 

In the terms, which couple more than two lattice variables, an analogous 
consideration leads to a 'maximal' factor 11^=1 cos Pi, where the momenta Pi 
may be all or partially different. The observation that variables can not be 
coupled over distances > 2n still holds. 

We have seen this restriction explicitly for n = 1. Since the couplings are 
confined to such a short range, we can easily include all of them. Unlike field 
theory, no truncation - which does harm to the improved properties of the 
action - is needed. This allows us to study the quality of a perturbatively 
improved action separately, whereas in field theory one can only study a 
superposition of the improvement and the truncation scheme. 

We also note that - within the limited set of couplings we deal with - 
locality becomes even better if the mass increases. For instance, a larger mass 
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Figure 4: The perfect 4-spin couplings oc X as functions of m. Couplings 
involving distances < 1 are on top, couplings involving distance 2 below. 

suppresses the couplings over distance 2 even more, relative to the leading 
coupling constants (see Table p. 

The restriction of the couplings to a finite range is in qualitative agree- 
ment with the effectively Id quark-gluon vertex function (quark fields con- 
stant in all but one direction) 0]. Moreover, also the increase of locality 
for rising mass agrees with fermionic models and with scalar fields in higher 



dimensions fB9f . 
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5 Numerical results for the energy gaps 



As a warming up exercise we computed the lattice partition function from 
direct integration. If we do so at A = and some A > 0, we can extract 
an estimation for -E'o(A) — Eq{0). In fact, this works much better for the 
perturbatively perfect action than for the standard action; for instance at 
m = A = 1 the continuum ground state energy Eq coui = 0.8038 is approxi- 
mated well using the improved action, 0.8021, whereas the standard action 
yields 0.7111. However, this involves additive constants, which may depend 
on A, so we focus on the energy gaps now. 

5.1 The simulation 

To compare the performance of the standard discretization and the perturba- 
tively perfect action, we simulated both actions at correlation lengths ranging 
from ^ f» 2 . . . 5, on a L = 30 lattice with periodic boundary conditions, using 
a standard Metropolis multi-hit algorithm. The first two energy gaps were 
extracted from the correlation functions (O|0(x)0(O)|O) and (O|0(a;)^0(O)^|O), 
and the statistical errors were estimated by jackknife analysis. The decay of 
the 0^ correlation function (in infinite volume) is given by 

oo 

(O|0(a;)V(O)^|O) = E |(O|0V)rexp(-AE„a;), (5.1) 

n=0 

(AE„ = En-Eo). For £ = 1, 2 it reduces to Q 

oo 

(O|(/)(x)0(O)|O) = 5]|(O|0|2n+l)pexp(-AE2n+ix), 

n=0 
oo 

(O|0(a:)V(O)'|O) = 5:i(O|0^|2n)pexp(-Ai?2nx), (5.2) 

n=0 

and for the harmonic oscillator we are left with 

(O|(/.(x)0(O)|O) = |(O|0|l)pexp(-AEix), 

(O|0(x)V(O)^|O) = KO|(/>^|O)p + |(O|0^|2)pexp(-AE2a:). (5.3) 

^At L < oo we actually obtain cosh functions, but we can easily measure the decay in 
a region, where this difference is negligible. 

^"^Due to the mirror symmetry of the potential, eigenfunctions for En have parity (— 
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In our simulations we study small interaction parameters A < 0.2. In this 
regime, varying the fitting ranges reveals that the i = 1,2 correlation func- 
tions do not pick up significant contributions from energy gaps higher than 
the leading ones given in eq. ( p.3|) . 

5.2 Numerical results 

Simulations were done for masses and anharmonic couplings in the range 
m = 0.2 . . . 0.5 and A = 0.001 . . . 0.2. We compared the first two energy gaps 
as "asymptotic scaling quantities", as well as AE2/AE1 as a scaling quantity. 
For a direct evaluation, we divide the lattice results by the corresponding 
continuum values. 

Figs. 1^ and ^ show the gaps AEi{X) and A£^2(A), measured at m = 0.5 
for the standard and perturbatively perfect action (and normalized by their 
corresponding continuum gaps AEn,cont)- These plots show that for anhar- 
monic couplings up to A ~ 0.05 the perturbatively perfect action reproduces 
the continuum gaps much better than the standard action. Comparing Figs. 
1^, P with Fig. ^ one recognizes roughly the same reliability range for first 
order perturbation theory in the continuum. This qualitative behavior of the 
improvement holds for a variety of masses. 

For the harmonic oscillator in infinite volume, the perfect action repro- 
duces the continuum gaps exactly, whereas the standard action yields 

AEn,stan{rn) = n ■ arccosh(l + m^/2). (5.4) 

(The value AEi^stanim = 0.5) ~ 0.495 agrees with Fig. |). The fact that 
even the standard action is perfect at A = is reflected by the exact values 
of the gap ratios. 

The scaling ratio ^§^(A) - again measured at m = 0.5 and normalized 
by the continuum value - is shown in Fig. |^. Unfortunately there is hardly 
any conclusive difference between the two types of action for that quantity. 

Nevertheless, in Fig. ^, which shows the gap ratios versus 1/^ at fixed 
A = 0.005, a better performance of the perturbatively perfect action is visible 
at intermediate correlation length. This is a scaling plot, based on 10^ Monte 
Carlo sweeps. 

As we know from Fig. the applicability of first order perturbation 
theory for AE2/AE1 is restricted to really tiny values of A. The hope that 
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Figure 5: The ratio AEi/ AEi^cont at m = 0.5, as a function of \. 



the perturbatively perfect action could still perform well beyond that range 
can not be confirmed. |^ In that range itself, the improvement is extremely 
difficult to demonstrate, because the standard action - being perfect at A = 
- is also excellent there. The available accuracy is mainly limited by the 
fitting precision of the exponential decays. 



6 Conclusions and outlook 

We have constructed a lattice action for the anharmonic oscillator, which is 
perfect to first order in perturbation theory. This is the first manifestly one 
loop perfect lattice action Comparing this action to the standard lattice 

^^We also looked at stronger interactions, such as A = 0.5, for masses m = 0.3 and 
m = 0.5. In that regime, the standard action scales even better than the perturbatively 
perfect action. 

^^For some time it was claimed that for asymptotically free theories, FPAs are not only 
classically perfect, but automatically also one loop quantum perfect. However, this claim 
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1.04 




Figure 6: The ratio l^E^j l^E^^cont at m = 0.5, as a function of \. 



formulation, we observe a clear improvement for the energy gaps AEi, AE2 
up to A = X/m^ ~ 0.2, over a variety of correlation lengths. As a conse- 
quence, (pseudo-) scaling laws of the type of eq. (|2.3| ) are improved at small 
interactions. In field theory, this corresponds to an improved asymptotic 
scaling. Unfortunately, for the scaling quantity AE2/AE1 an improvement 
could only be demonstrated laboriously. It seems to be restricted to very 
small values of A, in agreement with the performance of continuum pertur- 
bation theory. There the linear approximation is useful only for A < 0(10^^). 
One could have hoped that the improved action is successful also beyond this 
regime, but it turned out that this is not the case. This can be viewed as 
a negative sign for the direct application of perturbatively perfect actions, 
but the outcome might of course depend on the model. As a further test one 
could simulate the perturbatively classically perfect action for the 2d 0(3) 
model, which has been presented - but not tested - in Ref. |]13[ (Table 2). 

has recently been disproved [59. 
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Figure 7: The ratio ^^/^^^^,^^ at m = 0.5, as a function of X. 



Moreover, in the regime of tiny A, the standard action is also exceptionally 
successful in this toy model. Therefore, an accuracy of 4 or 5 digits is required 
to distinguish the results of the two actions in that regime. The fit of the 
exponential decay does not allow for such a high precision. The behavior of 
other scaling quantities, like AE^/AEi etc., is even worse, i.e. the regime of 
successful first order perturbation theory is even smaller. 

Actually this construction could be carried on to O(A^), but this includes 
couplings over distances 4, involving up to 8 lattice variables. Moreover, 
continuum perturbation theory suggests that the O(A^) perfect action would 
only help to proceed to slightly larger interactions, see Fig. 0. 



Still the anharmonic oscillator - and in particular the ratio of its mass 
gaps - may serve as a good testing ground for quasi-perfect actions, if we 
proceed to a non-perturbative improvement scheme Then one expects a 
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As a general trend, the perturbation series gets worse if higher gaps are involved. 
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Figure 8: The ratio AE2/AE1 at fixed A = 0.005. 



progress also for moderate and large interactions. Thus one overcomes the 
disadvantages of this model listed at the end of Section 1, although some 
truncation of the couplings will be needed. 

In particular, one may introduce an inverse temperature (3 in the ex- 
pression for the partition function plus RGT transformation term, and send 
/3 ^ 00. Then minimization is sufficient for a multigrid inverse blocking, 
and in this way one can identify a classically perfect action. Using the stan- 
dard action on the fine lattice and a fixed coarse configuration, we compared 
the minima for different blocking factors, and we typically observed a good 
convergence around blocking factor 10. As a test, we run the minimizer 
at small A and reproduced in this way the 4-spin couplings of Section 4. 

^"^For a finite blocking factor n one has to use the modified RGT parameter a„ = 
a(l - 

^'"^ Similarly, for the Schwinger model, the 0{g) (truncated) perfect plaquette couplings 
could be reproduce to percent level from the classically perfect action [|l^ (even though 
a slightly different RGT was used). 
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The 2-spin couplings bi arise from loop corrections, hence they are quantum 
effects, which are not present in the classically perfect action at small A. 

We can identify the classically perfect couplings over a wide range of 
interaction parameters. As an example. Fig. ^ compares Cq{X) for the per- 
turbatively perfect and for the classically perfect action. This figure does not 
imply that the perturbatively perfect approximation should be rather poor 
already at A ~ 0.1, as we have to conclude from the simulation results. 
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0.05 



0.5 



classically perfect 
perturbatively perfect 



1.5 



Figure 9: The dominant 4-spin coupling Cq at m = 1 as a function of the 
interaction parameter A for the perturbatively perfect action and for the clas- 
sically perfect action. 

Of course, the model is not asymptotically free, and the parameter A is 
relevant (not just "weakly relevant", i.e. in leading order marginal), so it 
does not belong to the class of models the classically perfect action is de- 
signed for. In fact, from the Table |I| we see that the bilinear contributions - 
that the classical approximation misses - are important. However, it might 
have a chance to perform well if we simulate at /3 > 1, where the quantum 
corrections to the perfect action are suppressed. The possibility that clas- 
sically perfect actions work to some extent also for models, which are not 
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asymptotically free, is conceivable and deserves being tested. 

Finally, thanks to the simplicity of the model, one can also perform the 
full path integral (numerically) instead, which yields a (quantum) perfect 
action. As a test, we also reproduced roughly some bilinear couplings 6j 
from Section 4 in this way. There one is restricted to small blocking factors 
and lattices, such that a good convergence of the action requires iteration. 
Then it is interesting to compare the couplings and their performance for the 
classically perfect action, and the action, which is - up to numerical errors 
and truncation - quantum perfect. The latter is more promising, but the 
classically perfect action is much of interest, because in complicated field 
theoretic models, this is what one - maximally - has at hand. 

An other important question is the convergence velocity in the multigrid 
procedure. In non-Abelian gauge theories, only very few iteration steps, with 
a small blocking factor (typically 2) are possible, hence a fast convergence 
to the FPA is crucial. We hope that starting from a perturbatively perfect 
action helps to accelerate the convergence. Also this can also be tested in 
the toy model discussed here. 

Furthermore one can use this model to test the "cycling" procedure of 
shifted forward and backward blocking, proposed by the Boulder group 



It is tractable in higher dimensions, but not strictly based on the renormal- 
ization group. Hence a toy model analysis of the errors emerging in the 
"cycling" process is of interest. 

At last - as we mentioned in the introduction - one may speculate that 
in complicated models the improvement could be pushed beyond classical 
perfection, by performing, say, one full block factor 2 RGT step at finite /3, 
starting from a classically perfect action. This method can also be tested for 
the anharmonic oscillator. 
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